A simple theory for high A/T c ratio in d-wave superconductors 
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We investigate a simple explanation for the high maximum gap to T c ratio found experimentally 
in high T c compounds. We ascribe this observation to the lowering of T c by boson scattering of 
electrons between parts of the Fermi surface with opposite sign for the order parameter. We study 
the simplest possible model within this picture. Our quantitative results show that we can account 
for experiment for a rather small value of the coupling constant, all the other ingredients of our 
model being already known to exist in these compounds. A striking implication of this theory is the 
fairly high value of the critical temperature in the absence of boson scattering. 
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A puzzling feature of high T c cuprate superconductors 
is the fairly high value of the maximum Ao of the su- 
perconducting gap compared to the critical temperature. 
Indeed it seems to range from 3 to 4 in most experiments, 
performed mainly on YBCO and on BSCCO @ . This 
is to be compared with standard BCS value 1.76 . Since 
it is widely believed that these compounds are uncon- 
ventional, with in particular changes of sign for the order 
parameter, it would seem that this is not much of a prob- 
lem. However this A /T c ratio is surprisingly very stable 
within all the generalizations of BCS theory which have 
been put forward for these compounds. Van Hove singu- 
larities and more generally any varying density of states 
raise it at most up to 2 , any reasonable anisotropy ||] 
gives a result not so much beyond the d-wave value 2.14 
and it requires strong coupling effects incompatible with 
experiments to push it in the experimental range. All 
these explanations are far from explaining the typical in- 
crease by a factor 2 compared to the BCS value, and one 
may wonder if a more complicated theoretical framework 
is not necessary in order to account for this ratio. 

We show in this paper that this is not the case and that 
the data can be explained quantitatively to a large ex- 
tent by a simple theoretical model in the standard frame- 
work of mean-field theory. Actually, except for the rather 
moderate value of the coupling constant we require, all 
the physical ingredients of our model are known to be 
present in these compounds, which makes our explana- 
tion a very natural one. However we do not apply this 
claim to the very high values of A /T c found recently 
in underdoped BSCCO. We believe that, even if they are 
clearly related to superconducting properties some 
additional physics, specific of this regime and this very 
anisotropic compound, is required to account for such 
very high results. Our focus is on the more standard val- 
ues found elsewhere, in other compounds (in particular 
YBCO) and in optimally and overdoped BSCCO. 

The basic idea of our model is the following. As is well 
known, when a superconductor has an order parameter 
which changes sign over the Fermi surface, superconduc- 



tivity tends to be destroyed by anything, like impurities, 
which scatters electrons between parts of the Fermi sur- 
face with opposite signs. If these scattering sources are 
present at T c but not at T = 0, they will lower the critical 
temperature but the zero temperature gap will be much 
less affected. This leads naturally to an increase of the 
Aq/T c ratio. In order for the number of these scattering 
sources to be temperature dependent, we have merely to 
take them as bosons, corresponding to a proper kind of 
collective modes of our system. Although other kind of 
fluctuations or modes may be considered, the simplest 
and most natural choice is phonon scattering. As we 
will see the typical energy needed for these bosons is in 
reasonable agreement with the frequencies available for 
phonons in these compounds. 

The idea of explaining a large value of Ao /T c by a de- 
crease of T c is already present in the literature, but to 
our knowledge it has not been put to work specifically 
in the case of d-wave superconductors. It is actually the 
usual qualitative picture behind the enhanced value of 
Ao/T c in strongly coupled standard superconductors : 
the argument is that thermally activated phonons tend 
to destroy superconductivity and lower T c while the zero 
temperature gap is not so affected since there are no real 
phonons present at T = 0. In the context of high T c su- 
perconductors Lee and Read Q, noticing the strong in- 
elastic scattering experimentally observed, have already 
proposed qualitatively this kind of mechanism to suggest 
a lowered T c . Here we rely specifically on the fact that the 
order parameter in high T c superconductors changes sign 
to obtain an important effect, compatible with experi- 
ment. We will more precisely assume the d-wave sym- 
metry, as it is most often done, although our mechanism 
actually requires only basically that the order parameter 
takes different signs on the Fermi surface. 

In order to explore this kind of explanation and see 
if it can work quantitatively for high T c compounds, we 
take the simplest possible model which retains all the 
qualitative features of our picture. Specifically we con- 
sider a class of models which has already been studied 
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by Preosti, Kim and Muzikar H in the presence of impu- 
rity scattering : we mimic a d-wave superconductor by 
taking an order parameter which takes a constant value 
A + on some parts of the Fermi surface and the opposite 
value A_ = — A + on the rest of the Fermi surface. We 
immediately specialize to the situation where the + and 
— regions have equal weight, as it is the case when they 
are related by symmetry. We assume that bosons scat- 
ter electrons from the + to the — region and vice- versa, 
and for simplicity we retain only those bosons. We take 
a simple Einstein spectrum with frequency Q for these 
bosons, with coupling constant A to the electrons. We 
assume the pairing interaction to have a characteristic 
energy much higher than T c and Q, and take a weak 
coupling description for the pairing. Therefore we do not 
make any specific assumption on the pairing mechanism : 
it may originate from pure Coulomb interaction or from 
spin fluctuations, or even have a more intricate origin. 
Again for maximum simplicity we keep only a constant 
repulsive pairing interaction between the + and the — 
regions. We do not expect any considerable quantitative 
changes from all these simplifications, all the more since 
it is known that the Aq/T c ratio is quite robust. 

With all these simplications the Eliashberg equations 
at temperature T read for our model : 
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Here A±^ n and Z±^ n are the order parameter and the 
renormalization function at the Matsubara frequency 
u> n = nT(2n + 1) in the ± regions. The effective 
frequency-dependent interaction A n = \ n — k contains 
the pairing interaction k , with a cut-off frequency oj c 
and the boson mediated interaction X p — Af2 2 /(S1 2 + cj 2 ) 
with Up — 2-KpT the boson Matsubara frequency. As 
mentionned above we have uj c ^> Q and T c . When we 
specialize to d-wave symmetry and insert the correspond- 
ing relation A_ „ = — A + n = A„ into these equations, 
we obtain Z_ .„ = Z+ <n = Z n and find that A n and Z n 
satisfy Eq.(l) and Eq.(2) (with A± ira and Z±, n replaced 
by A n and Z n ) provided that we change the sign in front 
of A„. The resulting equations are just the ones obtained 
in standard strong coupling theory, except that the roles 
are reversed between the boson and the Coulomb terms 
: the boson term is repulsive and the Coulomb one at- 
tractive. 

We have solved these equations directly both for the 
change of critical temperature and for the zero tempera- 
ture gap. However it turns out to be much more conve- 
nient to eliminate the pairing interaction and the cut-off 
in favor of the critical temperature T® in the absence of 
boson scattering. This is done by taking explicitely into 



account that, for fl <C lu 71 <C ui c , A„ goes to a constant 
and Z n goes to 1 . Let us first consider the calculation of 
T c , where A„ gets very small and Z n takes just its nor- 
mal state value. We call A^ this large frequency limit of 
A„ and set A n Z n = Aoo + d n . Since, as can be checked, 
the pairing term dominates in this range we obtain from 
Eq.(f ) and (2) (after taking into account that A„ and Z n 
are even functions of lu u ) that A^ satisfies : 



Aoo = 2/brT V 
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This leads to deal with S — 7rT^ d m /u> m Z m where the 
upper boundary can be taken as infinity since the sum 
converges. We have also to consider S' = nTJ2 ^/^mZ m 
where we have to keep the cut-off, but this can be ex- 
pressed in terms of T c ° as 5" = S z +(l/2) ln(T c °/T) + l/2/c 
with Sz = TrT^2(l/Z m — l)/uj m where again infin- 
ity can be taken as upper boundary. This leads to 
Aoo = S/[Sz + (l/2)ln(T c °/T)]. When this is carried 
into Eq.(f) this gives the numerically convenient eigen- 
value problem : 
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which is satisfied when T = T c . Note that a similar 
procedure could be applied to the standard Eliashberg 
equations. 

We are left with only two parameters, namely the re- 
duced frequency fl/T® and the coupling strength A of 
the bosons. We have plotted in Fig.l the ratio T c /T°, of 
the critical temperature T c compared to its value without 
coupling Xj?, for various values of the ratio O/T,? going 
from 0.2 to f . Naturally T c decreases with increasing A 
since boson scattering is pair breaking. High values of 
0/T° are not of much interest for us since they corre- 
spond to a full weak coupling regime, and the ratio A/T c 
will merely be given by the standard BCS value. Simi- 
larly large values of A lead to a strong decrease of T c as 
can be seen in Fig.f . This produces a large Vl/T c and 
leads again to the BCS value for A/T c . On the other 
hand the result in the low frequency limit £l/T® — s- is 
easily obtained. Indeed in this case uj n Z n = 7rT(2n+ f + 
A) and A„_ m + X n + m +i — A£ Iljm . This leads from Eq.(4) 
to an explicit expression for the n dependence of d n and 
to the result ln(T c °/T c ) = ip{\+ 1/2) - ip(l/2) where ip is 
the digamma function. This is just an Abrikosov-Gorkov 
result @, as might have been anticipated since bosons 
behave as impurities in the r2 — > limit Q. This ana- 
lytical result is actually in very good agreement with our 
numerical calculations for fl = 0.2. A noticeable feature 
of Fig.f is the strong sensitivity of T c on boson scattering 
even for moderate coupling strength. Indeed it is severely 
reduced already for A < 1. For example the slower de- 
crease of T c corresponds to the case il/T" — > 0, but even 
in this case the slope for small A is — 7r 2 /2. 
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FIG. 1. T c /T° as a function of the coupling strength A 
for n/T° = 0., 0.2, 0.3, 0.4, 0.6,0.8 and 1. Insert: A /T c ° as a 
function of A, for the same values of Q/T° ( = omitted). 

Let us consider now the calculation of the zero temper- 
ature gap, where the situation is somewhat more compli- 
cated. In the T — > limit A n and Z n become functions 
A(u>) and Z(u>) of the continuous variable u> n = lo. As 
for the calculation of T c we set A(oS)Z{lo) — Aoo + d(u>). 
In the same way the pairing term dominates in Eq.(l) for 
large u which leads to the following equation for Aoo : 
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where naturally A(u) in the right-hand side is expressed 
in terms of Z(ui) and d(uS) . For d(ui) we are left with : 
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where we can naturally use the even parity of A(u>) and 
Z{uj). As we have done above for T c we can use the 
weak coupling expression for the zero temperature gap 
Abcs = 1-76 T c ° for A = , obtained from Eq.(5) by 
setting A(uj) = A^ = AbcS: to eliminate the cut-off iv c 
and k. This leads to : 
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This last equation does not provide an explicit expression 
for Aoo m contrast with what we have for T c . However 
this is not in practice a problem, since it can be easily 
included in the iteration procedure used to solve numer- 
ically Eq.(6),(7) and (8). In order to find the gap Aq 



we still have to continue A(u>) and Z{ui) analytically to- 
ward the real frequency axis into A(y) = A{—iv) and 
Z{v) = Z(—iu), and solve A(A ) = A . This is done 
by using the explicit expression for this continuation [ pTj| . 
Note that this analytic continuation lowers noticeably the 
gap values, compared to the naive evaluation Ao = A(0). 

Just as for T c , the low boson frequency limit — > 
is of particular interest. Indeed it leads to a very sim- 
ple model of strongly interacting fermions with quite 
non trivial results. Naturally we have in this limit to 
let A increase in such a way that Afi stays finite oth- 
erwise one obtains trivially the BCS result, as it is 
clear from the equations found below. For this case 
Eq.(6) and (7) lead to algebraic equations because the 
lorentzian coming in the integrals gets very narrow. One 
obtains Z(w) = 1 + 7rAf7/2[w 2 + A 2 ^)] 1 / 2 and d{u) = 
-TTXnA(Lu)/2[u; 2 + A 2 (uj)} 1 / 2 , giving for A(u) the simple 
equation : 

AM 



A(u) = Aoo - ttAO 
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This equation (a fourth order equation) can be solved 
analytically and is very simple to solve numerically. The 
analytical continuation is merely obtained by replacing 
oj 2 by — v 2 in Eq.(9). However one can see easily that 
the equation A(i/)/Aoo = 1 - KA{v)/[-v 2 + A 2 ^)] 1 / 2 , 
where K = 7rAf2/Aoo, has no purely real solution when 
v/Aqo becomes larger than ^o/Aoo = (1 - if 2 / 3 ) 3 / 2 , 
corresponding to a gap A = A(i^o) = Aoo(l — K 2 / 3 ). 
Indeed, beyond this point, A{v) gets complex and the 
density of states is no longer zero. It is quite interesting 
to note that, in this limit, this nonzero density of states 
occurs before the equality A(Arj) = Ao is reached. To be 
complete we have to find the value of Aoo in this limiting 
situation. The integral in the defining Eq.(8) can actually 
be performed analytically (essentially by taking A(u>) as 
the variable) and is merely equal to ttK/A. Then A^ 
is obtained as the solution of the simple transcendental 
equation A^/ Abcs — exp(— 7r 2 Ai7/4Aoo) ■ Note that 
this equation has always a single solution in the physical 
range K < 1 . 

Coming back to our problem, we have, for fl/T® = 0.2 
and a varying coupling strength A, compared the re- 
sult obtained for Aoo in the limiting situation we have 
just considered with the general calculation we have per- 
formed for any Q,/T®. The agreement is quite good. 
From this it would be tempting to conclude that, for 
small we can also conveniently extract the gap it- 

self from this analytical solution. This is unfortunately 
not true : a good agreement for imaginary frequencies 
does not imply that the analytical continuations to the 
real frequency axis agree quite closely, since this analyti- 
cal continuation is very sensitive to small differences as it 
is well known. And indeed in our case there is a sizeable 
difference between the gaps obtained by the two methods. 
The results of our calculations for the zero temperature 
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gap are displayed in the insert of Fig.l. They show that, 
even for small A, A is quite sensitive to boson scatter- 
ing, although the effect is not as strong as for T c . We 
note in particular that the simple expectation that the 
gap would be essentially unchanged at T = because no 
bosons are present is not correct. 
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FIG. 2. Ratio of the gap over the critical temperature 
Ao/T c as a function of the coupling strength A, for fixed val- 
ues of n/T° = 0.2, 0.3, 0.4, 0.6, 0.8, 1. 

Finally our results for the ratio Aq/T c are shown in 
Fig. 2 for the interesting range of values for the parame- 
ter Q/T®. Wc note first that, for fl/T® = 1, the result 
does not depart much from the BCS result. Naturally 
this is even more so for higher values of 0/T° for which 
the results are not displayed. In the same way we find 
as expected that, for large values of A, A /T c decreases 
toward the BCS value. Because this is of little interest 
for our purpose, we have not explored further this regime 
where numerical calculations get more difficult since one 
has to deal with very different energy scales. The most in- 
teresting feature of our results is naturally the maximum 
obtained for A /T c at intermediate coupling strength. 
This maximum increases with decreasing O /T% while the 
A corresponding to the maximum increases at the same 
time. In particular for fl/T® = 0.2 we find A /T c close 
to 5. Clearly this trend continues as f2 goes to zero. In- 
deed as it is clear from above, the gap is independent of 
A in this limit whereas we have seen that T c decreases 
toward zero. Therefore we can in principle obtain a ratio 
Ao/T c as high as we like. However this would correspond 
to extreme parameter values. 

On the other hand we find among our results a range 
which is quite compatible with experiments. Indeed for 
£l/T® = 0.4 we obtain a broad maximum for A « .4 
with Aq/T c w 3.2. This fairly small value of A is quite 
reasonable. The value of A /T c is already quite con- 
sistent with experimental data, all the more if we take 
into account that anisotropy of the order parameter is 
likely to raise A /T c by itself, as it does for d-wave in 
weak coupling where this ratio is raised by 20 %. To 



be more specific, for fi/T c ° = 0.4, we find A /T c > 3.1 
for 0.36 < A < 0.52 . For a typical value of T c = 90 K 
, we find that the range of boson frequency goes from 
115 K to 170 K. For tt/T® = 0.3 we obtain correspond- 
ingly A /T c > 3.6 for 0.42 < A < 0.73, with Q rang- 
ing between 100 K to 170 K. This is a frequency range 
where an important weight for phonons is known to ex- 
ist in these compounds. Therefore, at least for optimally 
doped or overdoped compounds, our explanation for the 
high value of Ao/T c is completely coherent with exper- 
iment, which is quite satisfactory. For markedly under- 
doped compounds the general situation is no so clear and 
it is likely that the very high values observed in this case 
require an additional physical source which might for ex- 
ample be disorder. 

However a very striking feature of our interpretation 
is that it requires a fairly high value of T c °, that is the 
critical temperature without bosons, ranging from 290 K 
to 420 K for Vl/T® = 0.4 and from 330 K to 560 K for 
fi/T° = 0.3. Naturally it would be quite desirable to 
check experimentally this physical aspect of our model. 
One possible way would be to send a flux of phonons with 
the proper frequency to see if T c is affected as expected 
(these phonons could be generated themselves by tunnel 
junctions). Another much more interesting, though spec- 
ulative, possibility is to try to raise T c toward T® under 
static conditions. This could be done through a shift of 
the phonon spectrum, for example by applying high pres- 
sure in order to lower the number of phonons which par- 
ticipate in the decrease of T c . Actually it is known that, 
for Hg compounds, T c increases with pressure, which is 
compatible with our model. A test would be to measure 
also Ao under pressure and check that it is less sensitive 
to pressure than T c . This would be a clear indication 
that one can hope to raise T c even further by modifica- 
tion of the phonon spectrum. It is also important to keep 
in mind that it is often possible to obtain effectively an 
increase of pressure by proper chemical substitution in 
the compound. Therefore an understanding of the effect 
of pressure would open the way to a possible chemical 
increase of T c . Finally it is tempting to believe that the 
pseudogap observed above T c in undcrdoped compounds 
might be related to Tj? and could be obtained by treating 
our model beyond mean-field theory. 

In conclusion we have shown that the adverse effect on 
d-wave superconductors of boson scattering between re- 
gions with opposite sign of the order parameter provide 
a simple and natural explanation for the high values of 
A /T c observed experimentally. 
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